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Abstract 

The problem of quantum harmonic oscillator with " regular +random" square 
frequency, subjected to " regular +random external force, is considered in frame- 
work of representation of the wave function by complex- valued random pro- 
cess. Average transition probabilities are calculated. Stochastic density ma- 
trix method is developed, which is used for investigation of thermodinamical 
characteristics of the system, such as entropy and average energy. 

Introduction 

Recently there have been published a great amount of papers |1| concerned with the 
"quantum chaos", i.e. with the quantum analogues of classical systems possessing 
the dynamic chaos features. The investigations are conducted along the different 
directions, such as analysis of the energy levels distribution; definition and calcula- 
tion of quantities which are responsible in the quantum systems for the presence of 
chaos (corresponding to the classical Lyapunov exponents and KS-entropy); study 
of localization and delocalization of wave functions around classical orbits; etc. It 
is worth saying that though in most of the cases mentioned above one is faced with 



1 



the necessity to describe a quantum system statistically, so far there was not paid 
much attention to a stochastic behavior of the wave function itself. 

Many problems of great importance in the field of the nonrelativistic quantum 
mechanics, such as description of Lamb shift, spontaneous transitions in atoms, 
etc., remain unsolved due to the fact that the concept of physical vacuum has not 
been considered within the framework of the standard quantum mechanics. It is 
obvious that a quantum object immersed into the physical vacuum is an open system. 
There exist various approaches to the description of such systems, mainly in 
application to the problem of continuous measurements. One of them is based on 
the consideration of the wave function as a random process, for which a stochastic 
differential equation (SDE) is derived. But the equation is obtained by the method 
which is extremely difficult for application even in case of comparatively simple 
type of interaction between the system and the environment, so that some new 
ideas are needed Moreover sometimes it becomes necessary to consider the 

wave function as a random process even in closed systems (for example, when a 
classical analogue of the quantum system has the features of the dynamical chaos) 
i-i- 

To be able to describe the cases mentioned above, in the present paper we propose 
a radically new scheme of derivation of the evolution equation for a nonrelativistic 
quantum system, interacting with the thermostat (in particular with the physical 
vacuum), with the wave function represented by a complex- valued random process. 
The main idea of the new representation may be described as follows: a potential 
energy of the system "quantum object + thermostat" is assumed to be a random 
function. This is the case where the Schrodinger equation may be used only locally 
on small time intervals and may provide good phenomenological models for some 
problems, which are solvable within the framework of the multiparticle Schrodinger 
equation. 

In the majority of cases the main tool for investigating the particular problems 
is the perturbation theory, which sometimes fails to provide an adequate description 
of a real physical phenomenon. In the present paper the influence of the thermostat 
on an elementary process is considered nonperturbatively within the framework of 
randomly wandering one-dimensional quantum harmonic oscillator model. Average 
transition probabilities for the parametric oscillator are calculated, exact represen- 
tations are found for both widening and shift (analogous to the Lamb shift) of the 
energy level of the oscillator submerged into the thermostat (vacuum) as well as the 
entropy of an individual quantum state is calculated. 
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1 Description of the problem 

We shall consider the closed system "quantum object + thermostat" within the 
framework of a complex- valued probability process representation. The realization 
of the process is a wave functional '^stc{x,t\{(}), defined on L2 (^R^ ® "^{C}) ^P^<^® 
(extended space), where ((t) denotes a many-dimensional complex- valued random 
process. Time evolution of the wave functional is governed by the equation 

idt^^tc = H^stc, (1-1) 

where one-dimensional Hamilton function H is assumed to be quadratic over the 
space variable 

and the functions Q'^{t) and F{t) are random functions of time. Let them have the 
form 

F{t) = Fo{t) + yj2^2f2{t)Q{t - t2), (1.3) 

where ^l{t) and Foit) are regular (nonrandom) functions and /i(t),/2(t) are inde- 
pendent Gaussian random processes with the zero mean and 5— shaped correlators 

<Mt)f,{t')>=S,,6{t-f), I, J = 1,2. (1.4) 

Constants ej, i = 1,2 control the power of forces fi{t), t = 1,2. Functions Pi{t), i = 
1,2 are assumed nonnegative: pi,P2 > 0. Step-function Q{x) is defined by 

e{x) = 1 I ^ ^ ° (1.5) 

Let us assume that the following asymptotic conditions hold 

no{t) Fo{t) 0, Pi{t) 0, z = l,2. (1.6) 

t— >±oo V out ) t— >±oo t— >±oo 

which guarantee that the autonomous states 0^(x, t) exist as i ^ —00 

/ rr— \ 1/2 



is the wave function of a stationary oscillator and Hn{x) is the Hermi- 
tian polynomial. It also follows from ( |1.6|) that autonomous states 0°*'*(a;, t), which 



are obtained from ( |1.7] ) by a substitution of fij„ by Vtouti exist in the limit t -\-oo 
as well. B-functions in ( p..3|) reflect the fact that the random processes fi{t) and 
/2(t) are activated at the moments ti and ^2, respectively. If necessary, the functions 
Pi and p2 may be chosen having the form which prevents the jumps of VL and F when 
the noise is activating. The moments ti and ^2 are assumed to be finite to make 
the following inference correct. The aim of the paper is to find the average prob- 
abilities Wnm of transitions from the initial stationary states (f)^^{x,t) to the final 
ones t) when the evolution is governed by Hamilton function ( |1.2|) . Exact 

mathematical definition of Wnm will be given in the next section. 

2 Formal expressions for the wave functional and 
transition probabilities 



Proposition 2.1. The formal solution of the problem l \l.l\ )- ^T~^ ) may be written 
down explicitly for arbitrary fi^(t) and F{t). It has the following form 



^.tc(x,t|{C}) 



exp< I 



r](x - 7]) + —{x -riY + a 
2r 



X 



X — Tj 



r 



(2.8) 



where function x{y, t) satisfies the Schrddinger equation for a harmonic oscillator 
with the constant frequency Qin 

'dr 2d^y 2 ^' 

function rj{t) is a solution of the classical equation of motion for the oscillator with 
the frequency f2(t), subjected to the external force F{t) : 

7i + n\t)r] = F{t), r]{-oo) = f]{-oo) = 0, 

a{t) is a classical action, corresponding to the solution rjit) 



(2.9) 



a{t) 



n^rj^ + Fri 



dt\ 



[2.10) 



2 ' 2 

and r(t) and rit) are expressed in terms of the solution ^(t) of the homogeneous 
equation, corresponding to (\2.9{ ) 

'( + n^(t)^ = Q^ ^{t) ^ e^^-* (2.11) 
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as 

e(t)=r(t)e^W, rit) = m\, r = 7(t)Mn. 

There was introduced a special designation ( = for the set of functions ^(t) 

and ri{t) from ( ^.^ ). 

Proof. The proof is based on the substitutions first used in and may be per- 
formed by the exphcit verification under the only suggestion that all the executed 
manipulations are legal. A 

The set of solutions of type (2^) which is important for the following considera- 
tions in this paper is obtained from ( |2.8| ) after the substitution of xiu, t) by 4>T{y, t) 
from (|1.7|) . It is thus defined as 



^ite(a;,t|{C}) = -^exp|2 



f]ix - 1]) + —(x -riY + a 
Zr 

n = 1,2 



(2.12) 



The main properties of the set of functionals (|2.12|) , which are important in what 
follows are 

1. For any n the functional ^^"c (^j {C}) reduces to the autonomous state 0^"(a;, t) 
in the limit t ^ —oo. 

2. For any fixed (f elements of the set ( |2.12|) are mutually orthogonal in the sense 
of L2{R^), space of square-integrable functions: 

oo 

J ¥:2{x,t\{C})'^'kJ{x,t\{C}) dx = 6nm, (2.13) 

— oo 

where a bar denotes the complex conjugation procedure and Snm = 1, for n = m 
and Snm = for 72 7^ m. 

Definition 2.1. Average probabilities Wnm of transitions from the states \E'^"c^(x, t|{C}) 
to the stationary ones 0^*(x, t) in the limit t +oo are defined by 



*S(x,t|{C}) = E CnmmWZ\x,tl (2.14) 



m=0 

|2 



W„^= lim (|c„^r), (2.15) 
where the symbol < > denotes the procedure of averaging with respect to f\ and 

f2. 
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Definition 2.2. The generating function Istc{zi, Z2,t\{(}) for coefficients Cnm is 
defined by the expression 



Istc{zuZ2,t\{C}) = / dx'i'out{Zl:X,t)'^stciz2,X,t\{C}), 



(2.16) 



where 



oo 



SO that 



^stc{z,X,t\{C}) = E ^^i*i(x,t|{C}), 

oo n 

^out{z,x,t) = Y.^<P';:-'\x,t), 



Cnm(t\fl, f2 



dzfdz^ 



[2.17) 



21=22=0 



Proposition 2.2. Explicit expression for the generating function i\2.1d{ ) is given by 



IstciZi,Z2,t\{C}) 



1/2 



out 



\ ) 

where the following notations are made 
2 \ K 



exp{Azl+Bzl+CziZ2+Dzi+Lz2+M} (2.18) 



' 2Ue 



D = ^/2r]„„e^^°"** 



M 



^out [ ^out 



K 



1 



K 

iVto 



{QoutV -iv), 
if]' 



1] + 



K 



ut 



2K 



K 



K = -i^ + Qout- 

Proof. The proof is carried out by the direct summation of the series over Hermitian 
polynomials followed by the calculation of the Gaussian integral. A 
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Expanding the expression (|2.18|) in zi and Z2 powers, we obtain the coefficients 
Cnm- The ffist several of them are given below 

( '^y/^ifS^out \ M n T ^C^\ 

Coo — 777 e , Coi — i^Coo, cio — i^coo, l^-J-^J 



cu = {C + DL)coo, C20 = 72 (^5 + y ^ Coo, co2 = V2(^A + ^ 



Coo- 



Given formal expressions for the objects to be averaged it is necessary to reduce the 
averaging procedure to a form convenient for the subsequent analytical or numerical 
treatment. The following sections are devoted to the solution of this problem in 
different situations, which may occur when considering the general problem ( p..l|) - 

3 Average transition probabilities in case of ei = 

If ei = the function ^(t) from is nonrandom. We denote it as ^o(^), and get 

^o + nl{t)^o = 0, eo(t) ~ e^^-*, (3.20) 

eo(t)=ro(t)e^««=eoiW+^eo2(t). 
We also introduce the designation rjoit) for a function satisfying the equation 

iio + nl{t)7]o = Fo{t), r7o(-oo) = ?7o(-oo) = (3.21) 



Theorem 3.1. For any quantity G{'r]{t),ri{t)) local with respect to ri{t) and r]{t) 
(such are the coefficients Cnm in case of ti = 0) the averaging formula has the 
following form 



{G{'r]{t),r]{t))) = I I dxidx2 G{xi,X2)Pi{xi,X2,t\r]o{t2),'f]o{t2),t2), t > ^2- 

(3.22) 



where 



2vra. "^^\ 46163 -fei r ^ ^ 
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and 



hit) 



P2{t')eoi{t')dt' , hit) 



t2 



^2 



P2(t') (eoi(t'; 



t2 



262 



P2{t')Uit')U{t')dt\ 



t2 



yi = - 



io2{t)[xi -Vo{t)) -^02{t){x2 -r]o{t))] M„; 



Z/2 



Proof. It is obvious that the function 77 (t) is nonrandom in the time interval t < t2 '■ 
ri{t) = rjoit), while in the interval t > t2 it represents a random process with the 
evolution governed by the equation 



ri + n',{t)7] = Fo{t) + J2e2P2f2{t). 



(3.24) 



The initial condition for ( p.24| ) is defined by the requirement for the trajectory and its 
first derivative that they be continuous at the moment t2, i.e. 77(^2) = Voih), ^(^2) = 
7)0(^2)- To make the analysis of equation (|3.24|) , containing random processes, correct 
it is convenient to rewrite it as the set of two first order differential equations with 
respect to the quantities xi = r], X2 = rj 



Xi = X2, 

Fq - nlxi + V2e2P2/; 



X2 



2) 



Xi{t2) 
2:2(^2) 



^0(^2), 
^0(^2), 



(3.25) 



Equations (|3.25| ) are naturally interpreted as SDE for the random processes xi(t) 
and X2{t). Proceeding from them it is not difficult to write down the Fokker-Planck 
equation for the conditional probability density 



Pi(a;i,a;2,t|a;io, 0:20,^0) = {S{xi{t) - Xi)6{x2{t) - X2)) 



Xi(to) 

^2(^0) 



X20, 



describing the probability that trajectory {xi(t) , X2{t)) finds itself in the vicinity of 
the point (xi,X2) at the moment t, having started from the point (a:io,X2o) at the 
moment to- It can be shown that the equation for Pi has the form (see |]TD[ or ||11|| ) 



dPi 
~dt 



-X2 



dPi 
dxi 



(Fo - ^Ixi 



dPi 

dX2 



+ (^2P2 



d^Pl 

dxl 



(3.26) 
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The solution of (|3.26|) must be integrable and satisfy the obvious initial condition 



Pi = 5{xi - xio)5{x2- X20) (3.27) 

t=to 

As the position of the trajectory at the moment ^2 is known, it is natural to set 
to = t2, Xio = ?7o(^2)5 3^20 = Vo(t2)- The integrable solution of ( p.26| ), satisfying 



(|3.27| ), may be expressed in terms of the functions ^01(^)5 ^02(^)5 '7o(^)- It is easy to 
test that it has the form ( |3.23| ) by the explicit verification. The theorem is proved. 

A 

To find probabilities Wnm in this case we calculate the Gaussian integrals and 
then compute their limiting values at t ^ +00. 



4 Average transition probabilities in case of €2 = 

In the case of e2 = the time axis is broken into two parts. At t < ti functions r]{t) 
and ^{t) are nonrandom: rjit) = rjoit), C^{t) = C,o(t). At t > ti both r7(t)'s and ^(t)'s 
trajectories become random. 

Theorem 4.1. In the case of €2 = 0, at t > ti the set of equations l \2. !^ ), ( ^.lij) 



gives rise to the set of SDE, describing the evolution of four random processes 
u{t) = {u,{t),U2{t),U3{t),u,{t)) = (^7]{t),r]{t),Re (||) > (ft)) 

with the joint probability distribution function P2{u,t\uQ,ti), t > ti, satisfying the 
Fokker-Planck equation 

f) P 

= L2P2, (4.28) 

^ d 92 
L2 (m) = - ^ KiTr~ + ^iPi^i^ + ^iPi^ + 2eiPiUi- — h 4u3, 

^ OUi OU2 Oui OU2OU3 

Ki = U2, K2 = Fq- VlIui, Ks = ul - uj - -^4 = -2M3M4. 

and the initial condition 

4 

P2\t=ti = S {U - Uo) = YlSi^Ui - Uoi). 

1=1 

Proof. In the case under consideration the equation ( |2.9| ) is transformed to a set 
of SDE in the same way as it has been done in the previous section, namely, by 
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introducing the quantities Ui = rj, U2 = t) which reduce ( p.9| ) to ( p.25 ) with the only 
distinction in the initial condition: ■Ui(ti) = rjoiti), U2(ti) = r)o(ti)- 

The equation ( |2.11| ) may be reduced to a nonlinear first order differential equation 
by the substitution 



Uti)exp\'mt')dA, t>tu 



(4.29) 



which gives upon being applied to (|2.11|) the following SDE for $(t) in the interval 

tl < t < 00 



$(t) + $^(t) + n',{t) + J2eipifi = 0, ^tl) = Uti)/Uti 



(4.30) 



where the second equation expresses a condition which guarantees continuity of the 
function ^(t) and its first derivative at t = ti. The function $(t) is a complex- valued 
random process due to the initial condition. As a result the SDE ( [4. 301) is equivalent 
to a set of two SDE for real- valued random processes. Namely, introducing real and 
imaginary parts of $(t) 

$(t) = Usit) +iUi{t), 

we finally obtain the following set of SDE for the components of the random vector 
process u 



Ul = U2, 

U2 = Fq- njui - v/2eipin i/i, 

U3 = -uj + ul- Ql{t) - V2eiPi/i(t), 

U4 = — 2M3U4, 



Ml(ti) 
^2(^1) 

u^iti) 



Re {Uti)/Uti 
Im {Uti)/Uti 



(4.31) 



The pairs of random processes (ui,U2) and (usjU^) are not independent, because 
their evolution is influenced by the common random force fi{t). This means that 
the joint probability distribution 



Ui 



Uo 



P2iu,t\uo,ti) = {l[5iuiit) ^ 

U{ti) = Uo 

7]o{ti),r]o{ti), Re (Uti)/Uti)) , Im (Uti)/Uti. 



is a nonfactorable function. Proceeding from the known evolution equations (|4.31|), 



we obtain by the standard method the Fokker-Planck equation for P2 (see ||T0 
||1 1||), which has the form ( 4.28| ). The theorem is proved. A 



or 
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Given P2, one can average any quantity G {u{t)) which is local with respect to 
uit) : 



{G {u{t))) = J duP2 {u,t\uo,ti) G (u) , du = duidu2dusdu4. (4.32) 

But this formula fails to give a result if it is used for averaging the objects containing 
the coefficients Cnm from ( |2.17| ), which are nonlocal with respect to u. There does not 
exist a general approach to calculating the average value of any quantity nonlocal 
with respect to a random process. But it is known that for some types of such objects 
the averaging procedure may be reduced to finding a fundamental solution of some 
parabolic partial differential equation and its subsequent weighted integration. The 
description of the simplest case of this kind is given in |jl2|. It is not difficult to 
generalize the result obtained in |jl2|, and therefore the formulas (|A.4|) -( A.5 ) can be 
derived (see Appendix). Using ( |A.4| )- (|X75| ), we have the following proposition. 
Proposition 4.1. // the components of the random vector process u satisfy the set 



of SDE (\4.31\), then the averaging procedure can be represented as 



exp 



V,iu{T),u{t),t) dT~V2iuit)) 



due-^'^'^'^Q{u,u,t) , (4.33) 



where the function Q {u, u', t) is a solution of the problem 
dQ 



-^ = [L2 {u)-V,{u,u',t)}Q, 



Q (m, u', t) — > 5 {u — uq) , Q {u, u', t) 



(4.34) 



0, 



where 1 1 • 1 1 is a norm in R^. 

It is not difficult to show that when Vi = 0, the formula ([4.33|) transforms into 
( |4.32| ) with the substitution G = exp {— V2}, because in this case the equation ( [4.34| ) 



for Q transforms into the Fokker-Planck equation ( |4.28|) for P2. 

Using the proposition we obtain the representation for the average values 

^ +00. The exphcit 
1; m = 0, 1) : 



(|cnmP) , which are equal to the probabilities Wnm in the limit t 
expressions for the first four of them are presented here 



n 



0. 



nm 

{u, t) 



(4.35) 



where 

Hoo (m) 



out 



Mti)K\ 



exp 



^out ( 2Qout{u4^ + flout) 



- 2 
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(u) = 


UUrL I / 


Hio 


(u) = 




\Uti)K\' ^ 


Hn 


{u) = 







u% + 



\K\ 



-U1U2 



(4.36) 



K - {QoutUi - iu2){uiUi + i{u2 - U1U3)) Hoo {u) 



K = ilout + U4- ius, 
and functions Qnm (w, t) satisfy the equations 

9Qnn 



dt 



L'2Qnm ^nrriQr 



(4.37) 



Qnm (U, t) 6{U- Uo) , Qnm («, t) 

t >tl 



where 



Vnm = PnmU3, POO = POI = 1, PlO = Pll = 3 

To obtain Wnm it is necessary to proceed in ([4.35|) to the hmit t +00. The repre- 
sentation ( [4.35|) -( p:.37D is exact and free from any simphfying assumptions. Given a 
specific reahzation of pi{t), ( ^4.35| )-( |4737| ) is used as a basis for numerical calculations 
of the probabilities Wnm- 

More simple representation may be obtained under some additional assump- 
tions. First we introduce some useful designations. Namely, we shall start with the 
expression for the solution of equation (|3.20|) applicable at t +00 : 

6(t) = Cie*^°"** + C2e~^^°"'*, Ci = \Ci\e'^\ C2 = \C2\e'^\ 

then we shall write down also a corresponding representation for the solution of the 
equation (|3.21|) : 



Voit) 



i^od* + Cod), d{t) 



V^^in 



Ut')Fo{t')dt' . 



The following quantities, defined on the basis of the above formulas, will be included 
into the final expressions for the required probabilities 

C ^ 

' ci = lim d{t) = xA^e^^, (4.38) 



P 
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Theorem 4.2. Let the function Pi{t) have the form 

Piit) = 



1 ti <t <te, 
t>te. 



(4.39) 



and time te be large enough to guarantee that the solution Qnm {u, t) of the equation 
at the moment te may be replaced approximately by its stationary (att ^ +oo) 
limit. Then we can obtain the following representation of the probabilities Wnm 



where Qnm ^■^ solution of the shortened stationary equation 



L2 (O-Pnn^eS Q:L(e)=0 



(4.40) 



(4.41) 



Quantities Hnm for n = 0,1 and m = are given by 



- ' ^""*^^"[e3(ei+^i)-6-/i2]^^ 



exp 



-f^0l('Cl;6;'C3 



90 02 

Z,\Lout\L^n I , l2 ; 



2(6) 



'-[i2 - 66 - 6/^1 + /^2] i^oo(6, 6, 6 



2z/ 



fJ-i = -d^ + \ — (di cos P + d2 sin /3) 

V i ^i'n. 



2v 



S(6) 
+2 



in^ 'out 



= -d& + W — ((is cos (3 + di sin /3) 



{d{ + di^){l+ p) -2^ {di - di) cos 6 + 2rfirf2 sin 6 



(1-p) 

((i2C?4 + did^){l + p) + 2y/p {did4^ + d2d^) sin 5 + {did^ — ^2^4) cos 5 



(4.42) 

6+ 
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+ 



{dl + dl){l + p) +2y^ {dl- dj) cos 5 -2d3di sin S 



where 6 = 61 + 62- The ambiguity which arises in finding the solution Q^m 
equation (\i.4\ ) is eliminated by the requirement of the correspondence of Qnm 
stationary limit of the solution Qnm of the equation ( U-STj ). 

Proof. The proof originates from the representation ( 4.35|) -( 4.37 ). Denote the so- 
lution of the problem ( |4.37| ) in the time interval ti < t < t^ hj Q^m i^y ^) ^^^d that 
in the time interval t > t^ hy Qnm (^5 ^) • Then ("^5 ^) satisfies the first order 
differential equation 



dQ 



> 

nm 



dt 



1=1 



g VnmQlm, 



Qnm{u,t)-^Q<^{u,t,), Q>^(u,t) 0. (4.43) 



Using the method of characteristics the problem ([4 .431 ) can be solved for an arbitrary 
initial condition. So that we obtain the following expression containing the solutions 
of classical equations of motion ( |3.20| ) and ( |3.21| ): 



Qnm (U, t) = ( te) ■ CXp [ {A - p, 



U3 [^,t) dr 



(4.44) 



Quantities C,i{u,t) , i = 1, .., 4 in ([4.44] ) are the first integrals of the characteristic set 
of ordinary first order differential equations, corresponding to the equation ( |4.39| ). 
Their dependency on u variables is given implicitly by the relations 



U3 



■[6226 + 6216 + hi 



(eiiCa + ei2)(e2i6 + 622) + 611621^! 



U2 



■U4 



7:^[ei2^i + eii6 + h] 



(4.45) 



(6216 + 622)2 + 611^4 ' (6216 + 622)2 + 611^4 ■ 

The definitions of functions and constants appearing in (|4.45| ) are given by 

hit) = {d2de - c?44)6i(t) + (c/grfs - dide)^02it) + ^inVoit), ^2(i) = hi{t); 

621 (t) = di^02{t) - d2^0l{t), Cuit) = 622 (i), 

622 (t) = d4^oi{t) - d3^02(t), 611 (t) = 621 (t); 

dl = ^Ol{te), d2 = ^02(te), ^5 = r]o{te), 

ds = ioiite), di = ^02 (^e), 4 = VoiQ- 



(4.46) 
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The functions ^oi(^), ^02 (^) and rio{t) are defined in ( p.20|) and (|3.21| ). By the exphcit 



verification it may be proved that the new variables satisfy the initial conditions 

^i{u,Q = Ui, z = l,..,4, (4.47) 



so that the solution ( [4.44| ) satisfies the initial condition defined in ( [4.43|) 



Using ( [4. 451 ), the Jacobian 

of transition to new variables is calculated. In case of 1^4 = 0, which is important in 
what follows, it takes form 

= — ^ (4.49) 

C4=o (6216 + 622)^ ^ ^ 



Proceeding from the above we can rewrite the representation ( [4. 351) as: 



|Cnm|^) = J J Hnm{uii, t))Qnm (^5 ^e) exp ^(4 - Pnm) J ^3 (f, ^) dr 

(4.50) 

where integration over r is carried out under the assumption that ^ is constant. 
Assuming that the moment tg is big enough the solution Q^m {^j ^e) can be replaced 

approximately by the stationary limit Q^m ' which satisfies the equation 

(l^ - K^^) Q:1 = 0, = hm L, (4.51) 

Operator L2 is defined in ( [4.28|) . If necessary an exact value of tg, which will be 
sufficient for a legality of the replacement mentioned above, may be estimated. It is 
obvious that the stationary limit Qnmi^) function Qnm{u,t) depends on the 

initial condition for the latter. 

The insertion of the stationary solution Q^m into the formulas make them more 
simple, because the stationary equation allows a separation of variables. So substi- 
tuting 

Q':m = m)Qnmi^l,^2,^3), (4.52) 

into ( [4.51| ), we obtain the following (shortened) equation determining the function 

- PnUs) QLiO = 0, (4.53) 
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out y 



Substituting ( [4.52| ) into ([4.50|) and using (^4.49|) , we come to the formula 

= f^rr / rfeirfe2cie37%%^^^n^ (^(e, t)) i^.^o- (4-54) 

J 691 (ta + 699 F"™ ^ / ?4 ^ 



To find the average probabihties Wnm we proceed to the hmit t +00 in ( |4.54| ). 
This may be done for any n and m. After implementing simple but wearisome 
transformations the limiting values of ( [4 .541 ) were obtained for the minimal n, m. 



They have the form ( |4.40| )- (|4.42| ). The theorem is proved. A 



It is a real challenge to obtain the solution Q^mi corresponding to a given initial 
condition appearing in ( |4.37| ), without solving the equation on the whole interval 
ti < t < 00. But this problem may be resolved in the following particular case. 
Theorem 4.3. Let Fo{t) = const = Fq, Qo(t) = const = Qq, then Q^m entering 
into the statement of the theorem is given by 

QnmiUl, M2, U3) = U2, U3), (4.55) 

where g^*^(Mi, M2, M3) is an arbitrary fixed solution of ( H-SSj ) decreasing on infinity, 
and the constant Cnm is expressed in terms of this solution as follows 

Cum = >;t(Mo) |y duidu2du3 («) [^^^g^L («! , ^2, ^s) } (4.56) 
where is any solution of the conjugate equation 

[^'Y (u) - Pnmus) Y:Uu) = 0, (4.57) 



- st\ + 4 Q (p- (p (p 

Kt' = m_, A'f = -fiL"i. A? = "5 - - f!L), A'f = -2113114. 

Proof. To prove the above statement we consider the quantity 

lZ!{t) = Jdu Ynmiu, t)Qnm{u, t), (4.58) 
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where Qnm{u, t) is a solution of the problem ( |4.37|) and function Ynm{u, t) is assumed 



to be found proceeding from the requirement that ([4.58|) does not depend on the 
time variable, i.e. dl^/dt = 0. Integrating by parts we can easily show that the 
requirement will be satisfied if the function Ynm is an arbitrary solution of the 
equation conjugate to ( [4.37| ) 



^ = {-Lt + Vnm) (4.59) 

Lt (u) - E^^l" + ^1^1-?!^ + ^1^1 III + 

The only condition which is supposed to be imposed on Ynm is its sufficiently rapid 
decreasing with ||-u|| —>■ cxd. This guarantees that all the terms, containing integra- 
tion over infinitely remote surfaces, will vanish. 

It is worth mentioning that the quantity (|4.58 ) does not depend on time even 



when the functions ^o{t) and -fo(^) entering into the statement of the theorem are 
nonconstant. This condition becomes important in the case when the constancy 
of the quantity (|4.58|) is used for the restoration of information about the initial 



condition, which is lost in proceeding to the stationary limit in ( [4.37| ). 

If pi(t) = P2it) = 1, -fo(t) = const = Fq, Qo(t) = const = flo the operator 
which is conjugate to the Fokker-Planck operator L2 does not depend on time. In 
this case it may be the time independent solution Y^^ of the equation ( [4.59| ) 

(-L++p„„M3)r„t = 0. (4.60) 



to be substituted into ( [4.58| ) instead of Ynm- Note that in this case L2 = (-^2*) 



so that the equation ( [4.60|) is in fact a stationary equation corresponding to ( |4.59| ). 
Having chosen such a function Ynm one can equate values of the quantity ([4.58|) at 
t = ti and at t ^ 00 : 

i;t(no) = / duY:Uu)Q:Uu). (4.61) 

Proceeding from an arbitrary solution Qnmi^) '^^ the stationary equation (|4.51| ) with 
a sufficiently rapid decreasing on infinity we obtain the required solution Qnmi^) 

Q'^u) = CnmqLiu), (4.62) 

where the constant Cnm depends on the specific choice of q^^{u). Unfortunately, so 
far we do not possess a strict proof of such representation, but it seems plausible in 
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view of the results obtained in exploring the more simple one-dimensional equation. 
Substituting (|]6|) into (^381) , one readily obtains 



Taking into account that variables in stationary equations are separated, i.e. Q^j^{u) - 
5{uA)Qnm{'^i, U2, U3) q^mi'^) = <5(M4)g^m(^*b ^2, U3) , the formula (14.561) for Cnm is ob- 
tained. A 



5 Average transition probabilities in case of ei 7^ 
and 62 7^ 

If influence of both random forces is taken into account, i.e. ei 7^ and 62 7^ 0, 
evolution of the system depends on the correlation of the moments ti and ^2- If 
t < ty = max{ti,t2) equations ( p.9|) and (|2.11|) for the trajectories ri{t) and ^(t) are 



different, depending on whether i> = ti or t> = t2- The two possibilities correspond 
to the two probability distributions of random variables T]{t^) and ^(t>), which are 
the initial values for the corresponding trajectories on the interval t > t>. Denote 
Zi =77, Z2 = f], Z3 = Re{^/^), Z4 = Im{^/^). For the probability density function 
R{z,t), which is defined as a probability for the trajectory z{t) to be found in the 
interval [z, z + dz\ at the moment t > = min(ti,t2), we can write down the 
following expression 

{Pi {Zi, Z2, t-^\zoi, Z02, t<)S{z3 — 2:03)^(2:4 — ^04), t> — ti, 
(5.64) 
P2{^,ty\zQ,t^), = t2, 

where 

zo = (^01,^02,^03,2:04) = {vo{t<), ?7o(t<), Re{io{t<)/^o{t<)), /m(^o(t<)/6(^< 



and functions Pi and P2 are the solutions of (|3.26|) and (|4.28| ), respectively. It is 



obvious that the normalization condition / R(z,t^)dz = 1 holds. 

At t > t> we have the same set of SDE for the components of the random vector 
process z{t) both at t> = ti and t> = ^2- Its inference literally reproduces the 
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derivation of the set (^4.31 ) and results in 



Zl — Z2, 

Z2 = Fq- Vtlzi - V2eipi2;i/i + ^/2e^f2, 

= Zl - Zl - ^l(t) - v/2^/l(t), 



(5.65) 



z{U 



z>, 



where distribution R{z-^,ty) of the components of the random vector i> is given by 
the formula (|5.64| ). The representation of the joint probability density 



P3{z,t\zy,ty) = {6iz{t) - ^)) 



derived by the standard method from (|5.65|) , is given by the Fokker-Planck equation. 
Thus we arrive at the analogue of the theorem 

Theorem 5.1. In the case of ei ^ and 62 ^ the set of equations \2. 9[ j, i j| j 
generates att > t-^ the set of SDE describing the evolution of four random processes 



z{t) = {z^{t),Z2{t),z^{t),z^{t)) = {riit),f,{t),Re 



At), 



Im 



with the joint probability distribution Ps{z,t\z^,t^), t > t^ satisfying the Fokker- 
Planck equation 

' LsPs, (5.66) 



dt 



Uiz) 



^ d 



92 



92 



i=i 9z, 
Ki = Z2, K2 = Fo 
and the initial condition 



+ {e2P2 + ^iPizt)-Q^ + ^^P^'g^ + 2eiPi2i 



Za Z-i f^Q, 



dz2dz3 

- —2ZsZ4. 



+ 4^3 



p. 



6{z- z>] 



with the probability distribution of components of the random vector i> given by the 
formula l \5.64\ )- 

It is not difficult to show that if ei = the equation ( |5.66D transforms into ( |3.26D 
by the substitution 



P3{z,t\z>,ty) = 6 [z3- Re 



Ut)/ 
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5 \ Z4 — Im 



P{zi, Z2,t). 



When €2 = the equation ( p.66| ) transforms directly into ( 4.28 ). Thus the relation of 
the general case with the particular situations, considered in the previous sections, 
is established. 



Using again the formulas (|A.4|) - (|A.5|) , the representation (|4.35| ) of (|c„mP) may 



be generalized to the case under consideration 

'^^ ^ J dz>R{z>,ty) j dzHnm{,z)Qnm{,z,i). (5.67) 



nm 



As compared with ( ^.35| ), an additional integration with the weighting function 



R{zy,ty) allowing for the dispersion of initial values of the trajectory, was included 
into ( |5.67] ). Functions Hnm{z) are defined in ( [4. 361 ) and functions Qnm{z,t), depend- 
ing on z^ as on a parameter, are solutions of the following problem 

^Qnm _ f ^ _\/n 
Ot 

Qnm{z,t) 6{Z- Z>), Qnm{z,t) 0, (5.68) 

where 

Vnra = PnmZs, PoO = Poi = 1, Pw = Pll = 3. 

Proceeding from (|5.67|) - (|5.68|) the theorems analogous to the theorems and ^.|^ 



can be proved. Here we shall combine them into the one 

Theorem 5.2. Let floit) = const = Qq, Foit) = const = Fq. Let the functions Pi(t) 
have the form 

and let the moment te lie far enough in the future to provide the legitimacy of the ap- 
proximate replacement of the solution Qnm {u, t) of the equation ( H-STj ) taken at this 
moment by its stationary limit att ^ oo. Then we have the following representation 
for the probabilities Wnm 

Wnm = n\rA^m{t>) j d^.d^^d^.f^U^,, ^2, ^3)Hnm (6, 6, 6) ■ (5.70) 
^nmit^) = J dz^R(z^,t^)Cnmiz^) 

where Hnm are defined in l\4-43i )) the function q^^ is an arbitrary (decreasing on 
infinity) solution of the shortened stationary equation 

h -PnU3)qL = 0, (5.71) 
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- - E + (^2 + .l^S^ + J + 2e,?,g|^ + 2.3. 

and t/ie function Cnm{z>) depending on the parameter z> is defined by the formula 



Cnm — Ynm{^>) dzidz2dz3 Y^^{z) ^^^^q^^{zi , Z2, Z3) j (5.72) 
where Y^^{z) is an arbitrary solution of the conjugate equation 

(^{LfY (z) - pr^mzs) Y:'„X^) = 0, (5.73) 

irf = Z2, irr = -fiL^i, K^,' = zl-zl-nlJ, Kt = -2z,z,. 



6 Another approach to the calculation of average 
transition probabilities 

There is an alternative approach to calculation of the average transition probabilities. 
The main idea of the method will be demonstrated with the particular case of the 
problem (|L1]) - (^. 
Theorem 6.1. Let in ( \TJi ) fio = 

= 0, ei = 0. Then the solution of the problem 

j| j-( |7!^ may be represented as 

^stcit, x) = e^(«i(^-«2)+«3)^(^ _ (6 74) 

where function (j){y, t) satisfies a usual partial differential equation 

d(j) 1 d'^cj) i 2 
'Tt=-2-¥-y^2y 

and random processes ^i, ^2, ^3 satisfy the system of ordinary SDE 



dii{t) = -i^2it)dt + ^/2^dW2{t) 
d^t) = Ut)dt 

dUt) = l{ei{t)-ze{t))dt + V2^2Ut)dW2{t). 
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(6.75) 



Proof. Let us rewrite the equation (|1.1| ) clS cL SDE 



11 



192^ 



stc 



stc 



2 dx^ 



dt - y/2^e{t - t2)x'^stcdW2it), 



(6.76) 



which determines an increment of the random process "^stc during the time interval 
dt. The quantity W2{t) in ( |6T6|) is the Wiener process (W2{t) = Jl f2{t')dt'^ , i.e. 
the Gaussian random process which is completely determined by the initial condition 
^^2(^2) = and transition probability 



P{wo,to\w,t) 



1 



27l{t - to 



: exp 



2{t - to) 



The equation ( |6.76| ) differs from an ordinary SDE in that it includes derivatives of 
the process with respect to parameter x. The idea of the method is to represent 
\E'stc as a function of several random processes satisfying the set of ordinary SDE, 
the latter being chosen so as ( |6.7(j| ) would be turned into an identity after the 
substitution of "^stc- In doing so it is necessary to apply the Ito calculus as follows. 
If "if stc = 0(0; "where ^ = (^1, ...,^n) is a set of random processes satisfying the set 
of ordinary SDE 



dUt) = aiil t)dt + h{l t)dW2{t), I = 1, n, 



(6.77) 



and is a differentiable function of its arguments, then an increment of the process 
^stc is 



d'ifstcit) 



+ T.b^J^dW2{t). 

i=l '^'^i 



(6.78) 



where d^i are infinitely small increments of the processes determined from ( |6.77| ). 
Using this rule, we can show after simple but cumbersome manipulations that the 
equation ( |6.76| ) holds at t > ^2, if "^stc is represented by ( |6.74| )-( |6r75| ). A 

The initial conditions for may be set to zero: ^1(^2) = ^2(^2) = ^3(^2) = 0. 



Then solving the equation 



at t < ^2 and denoting the solution as "^ftdt, x), we 



obtain the initial condition for as: (f){t2,x) = '^f^^{t2,x). The boundary conditions 
for are imposed by the natural requirement of normability of the solution. 
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Having the solution ( 6.74 ) it is not difficult to find the matrix elements of dif- 
ferent quantum mechanical operators and the probabilities of transitions into given 
stationary states (for example, to (p^^)- To be able to average these quantities over 
the joint distribution P{t, C,), "we must solve the Fokker-Planck equation correspond- 
ing to the set ( |6.75|) . 



Unfortunately a substitution reducing the equation (|1.1| ) to the set of ordinary 
SDE is unknown in the general case (|L1])-(|L6|). Therefore at present the described 
method provides noting more than an instructive illustration of the one of possible 
approaches to the problem to be solved. Though not leading to the final results, it 
reserves the field for the further study. Undoubtedly it would be very interesting to 
obtain the solution with the help of the described procedure, even if only to compare 
it to those found in the previous sections. 



7 Thermodynamics within the framework of rep- 
resentation by stochastic density matrix. Ther- 
modynamical characteristics of oscillator 



It is well known [|T3| that the key object of interest in quantum mechanics is the 
density matrix. 

Definition 7.1. The stochastic density matrix is defined by the expression 

oo 

Pstc (x, t; {C}|x', t'; {('}) = Y: ^'"Vi™^ {x, t; {C}\x', t'; {('}) , (7.79) 

m=0 



pi^) {C}\x'A {('}) = J^^irc^ (x,t|{C}) vl/i^ {x',t'\{C}), (7.80) 

where has the meaning of the initial distribution over quantum states with ener- 
gies = + ^) ^ini until the moment when the generator of random excitations 
is activated. 

Definition 7.2. The expected value of the operator A {x,t\{C,}^ in quantum state 
with the index m is 



(7.81) 
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The mean value of the operator A (^x,t\{(}^ over the whole ensemble of states will 
respectively be given by 



A = ^hm^Sp:, Sp^^^Apstc I Sp^ Sp^^^pstc } 



(7.82) 



The operation Spx in (|7.81|) and (|7.82|) is defined by 



Spx{K{x,x')} 



TT 



dxK{x, x) 



(7.83) 



for any function K{x,x'). 

Using (|7.80|) and the properties of the functionals \l/^tc^ {x, ^KC}) "we easily obtain 
the expression for the total nonstationary distribution function 



Wo 



0<J 

SpxSp^^y {pstc {x,t; {C}\x',t'; {('})} = 



(m) 



(7.84) 



m=0 



If the initial weighting functions Wq™^ given by the canonical distribution 



(m) 



exp {—Em/^T), the expression ( [7.84] ) takes the form of the Planck distribu- 



tion (see [0) 



Wo iP) 



f3 



(7.85) 



e/3 _ 1 ' - kT ■ 

Substituting the expansion (|2.14| ) of the wave functional in out-states into ( [7.79|) - 
(|7.8CI| ) we have the following representation 



Pstc {x,t;{C}\x',t';{C'} 



E ^t^^^k {t\{c}) c^i {t\{Q)<PT\^.t)<pr\^'.t) 



(7.86) 



m,k,l=0 



Definition 7.3. The nonequilibrium partial distribution function is defined by 



w 



(m) 



ei,e2,t) = Spr^y \ E 



.k=0 



(k) 
Wo 



(^kr, 



- w. 



(m) 



Cmk [t\{C,} 



I (™) 
+ ^0 



CXJ 

E { A,„ (t) - 4"^) (t)] } + 4' 



fc=0 



(m) 



(7.87) 



where 



^km {t) = Sp^^^ Ckm {t\{C}) = ( Cram {t\{Q) 
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In this case the total distribution function is equal to the sum 

oo 

Wo= E^^"Uei,e2,t). (7.88) 

m=0 

In case under consideration one can introduce different definitions for such ther- 
modynamical quantity as an entropy. Despite formal similarity definitions done 
may provide or not the connection of defined quantity with irreversibility of the 
system evolution. For example one can define the total and the partial entropy in 
the following way. 

Definition 7.4. The formal total entropy of nonequilibrium state is defined as 

S (ei, 62, t) = -Sp^^ySp^ {pstc^T^ Pstc} ■ (7.89) 



Definition 7.5. The formal partial nonequilibrium entropy is defined as 

(6i,62,t) = -Sp^^^^Sp. {pi^ Inpi^ } . (7.90) 

It is not difficult to show that the formal partial entropy does not depend on time 
and has no relation to thermodynamical irreversibility. 

Proposition 7.1. For any m the formal partial entropy S^"^^ (ei,e2,t) is equal to 
zero. 

Proof. Let's consider the A^-dimensional square matrix A with elements = aiflfc, 
where a^, i = 1, . . . , are the elements of iV-dimensional vector. It is possible to 
find all eigenvalues Aj and to find out the structure of eigen-subspaces for matrix 
A. Namely, one can show that Ai — a\ -\- a\ -\- . . . -\- a%, X2 — X3 — ■ ■ ■ — Xn — 0. 
At that eigenvector ei coincides with a, and eigen-subspace corresponding to zero 
eigenvalues is orthogonal to a. 

Generalizing this result on the case of infinitely dimensional matrix p^^J , one ob- 
tains: there is one eigenvector "^^^^J (^x,t\{(}^ , corresponding to nonzero 

eigenvalue Ai = Spj^^^Sp^ \^p^stc} 7 ^-iid there is an infinitely dimensional eigen- 
subspace, corresponding to zero eigenvalue, which is orthogonal to this vector. Sup- 
plementing the vector (7r/Jli„)^^^ "^^J^J (x, i|{C}) with any orthonormal set of vectors 
lying in the subspace mentioned above, one obtains the basis of the whole space 
which brings matrix p^JJ^ to diagonal form. Understanding the uncertainty In as 
a limit 

OlnO = limslns = 0, 
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one obtains for formal partial entropy: 

Sf^^ (ei, £2, t) = Sp^^Sp^ {piS} ■ Sp^^ySp^ {inpSS } = 0, 

which makes the proof complete. A 

If one wishes to have the quantity describing irreversible behavior of the system, 
it is necessary to change definition of entropy. 

Definition 7.6. Total and partial entropies of nonequilibrium state are defined as 

S {€i,€2,t) = -Sp^ {Pav^^Pav} ■ (7.91) 

and 

(61, 62, t) = -Sp, {pi-) lnp(-)} (7.92) 

correspondingly, where 

Pav = Sp^^y {pstc} , pIT^ = Sp^^y {p'^^J } . 

Unfortunately we have no at the moment simple enough analitical representation 
for the quantities defined in such a way. 

To illustrate the definitions given above we calculate the average energy of oscil- 
lator in the ground, vacuum, state (i.e. at m = 0) assuming that both regular and 
stochastic parts of the external force are absent. In this case the density matrix has 
the form 

pTc (x, t- {Q I x\ t'- {C'}) = exp [-^{x' + x'^)- 

> 

(7.93) 



1 I- 



2 



J $ (r) d^-\j ^* (^) dr-\[^ it) - {t') ^ 



Proposition 7.2. Let 0(i) = Qi^, -Fo(^) = 0, 62 = 0, P\{t) = 1. Then the average 
energy 

Eif^W = Sp^Sp^^y (^pS) 

is represented by 

41 (A) = 2 1 - 7v / dusK, (A, «3) + ^ / dusK2 (A, %) [ , (7.94) 



with the designations 



Ki (A, Us) = CoMs?"* (A, Us) X 



26 



^0-1 



Us 



2^00 



00 L 



' ^00 + 1 ^ Us I Aqo 

2^00 



00 



'A 



00 



K2 (A, Us) = Cousq'^ (A, Us) X 
Us /a 



^^00 



00 



00 



1 _^ M3 / Aqo 



2A§o 



-^^00 



Aqo (A, U3) 



A 



The function q^^{us) is an arbitrary solution of the equation 



d 1 / 2 



decreasing as \us\ 



^3 dus 
00, and the constant Co given by 

Y'\uso,U4o) 



Co 



Jdus Y't{us,0)q''{us)' 



Here the function Y'^^{us, U4) is an arbitrary solution of the equation 



2vst 



d'Y 



ei- 



dul 



- [Un 



Ua 



^yst 
dus 



dY^^ 

2usU4- UsY'^ 

dUd 



00. 



(7.95) 



(7.96) 



(7.97) 



(7.9^ 



decreasing as u| + ul 
Proof. In fact the proof copies the manipulations performed repeatedly in this paper 
and thus may be omitted. A 

The second term inside the figure brackets in ( [7.94[ ) is a level shift which is well 
known from quantum electrodynamics as the Lamb shift, the third term determines 
the magnitude of the ground state energy broadening. Note that the lifetime at this 
level is proportional to the inverse of the broadening 



H-cxo 



dusK2 (A, Us) 



(7.99) 



The average energy of a randomly wandering (QHO) for any quantum level is 
calculated quite similarly. 
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Conclusion 



There are three different reasons which may cause a chaos in the basic quantum 
mechanical object, i.e. the wave function. The first reason refers to measurements 
performed over a quantum system The second reason consists in the more 

fundamental openness of any quantum system resulting from the fact that all the 



beings are immersed into a physical vacuum [T^]. In the third place, as it follows 
from the recent papers P, |3, ||, a chaos may also appear in the wave function even 
in a closed dynamical system. As it is shown in 0, there is a close connection be- 
tween a classical nonintegrability and a chaos in the corresponding quantum system. 
Many of the fundamentally important questions of the quantum physics such as the 
Lamb shift of energy levels, spontaneous transitions between the atom levels, quan- 
tum Zeno effect |[T^, processes of chaos and self-organization in quantum systems, 
especially those where the phenomena of phase transitions type may occur, can be 
described qualitatively and quantitatively in a rigorous way only within the nonper- 



turbative approaches. The Lindblad representation |T5|, for the density matrix 
of the system "quantum object + thermostat" describes a priori the most general 
situation which may appear in the nonrelativistic quantum mechanics. Nevertheless, 
we need to consider a reduced density matrix on a semi-group 0], when investigat- 
ing a quantum subsystem. This is quite an ambiguous procedure and moreover its 
technical realization is possible only in the framework of a particular perturbative 
scheme. 

A crucially new approach to constructing the quantum mechanics of the closed 
nonrelativistic system "quantum object + thermostat" has been developed recently 
by the authors of 0, |^ from the principle of "local correctness of Schrodinger rep- 
resentation". To put it differently, it has been assumed that the evolution of the 
quantum system is such that it may be described by the Schrodinger equation on 
any small time interval, while the motion as a whole is described by a SDE for the 
wave function. In this case, however, there emerges not a simple problem to find a 
measure for calculating the average values of the physical system parameters. Nev- 
ertheless, there exists a certain class of models for which all the derivations can be 
made not applying the perturbation theory 0. 

In the present paper we explore further the possibility of building the nonrela- 
tivistic quantum mechanics of closed system "quantum object + thermostat" within 
the framework of the model of one-dimensional randomly wandering QHO (with a 
random frequency and subjected to a random external force). Mathematically the 
problem is formulated in terms of SDE for a complex-valued probability process de- 
fined on the extended space R^^R^^y. The initial SDE is reduced to the Schrodinger 
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equation for an autonomous oscillator defined on a random space-time continuum, 
with the use of a nonlinear transformation and one-dimensional etalon nonlinear 
equation of the Langevin type defined on the functional space R{^}- It is possible to 
find for any fixed {^} an orthonormal basis of complex- valued random functionals 
in the space L2{R^) of square-integrable functions. With the assumption that the 
random force generator is described by a white noise correlator, the Fokker-Planck 
equation for a conditional probability is found. From the solutions of this equation 
on an infinitely small time interval a measure of the functional space R^^y can be 
constructed. Then by averaging an instantaneous value of the transition probability 
over the space R{q, the mean value of the transition probability is represented by 
a functional integral. Using the generalized Feynman-Kac theorem, it is possible 
to reduce the functional integral in the most general case, where both frequency 
and force are random, to a multiple integral of the fundamental solution of some 
parabolic partial differential equation. The qualitative analysis of the parabolic 
equation shows that it may have discontinuous solutions This is equivalent to 
the existence of phenomena like the phase transitions in the microscopic transition 
probabilities. In the context of the developed approach the representation of the 
stochastic density matrix is introduced, which allows to build a closed scheme for 
both nonequilibrium and equilibrium thermodynamics. The analytic formulas for 
the ground energy level broadening and shift are obtained, as well as for the entropy 
of the ground quantum state. 

The further development of the considered formalism in application to exactly 
solvable many- dimensional models may essentially extend our understanding of the 
quantum world and lead us to the new nontrivial discoveries. 

Appendix 



Theorem A.l. Let the set of random processes {C,i,C,2, ■■,^n) = ^ satisfy the set of 
SDE 

n 

ii{t) = ai (f, t) + XI ^ij fj{t), i = 1, n, 

(/.W/.(0) = M(i-0, 

so that the Fokker-Planck equation for the conditional transition probability density 

p^'^ (6,t2ifi,ti) = {s {m) - 6)) t2 > ti (A.l) 
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is given by 

ap(2) _ 
at ~ 



(A.2) 



The processes are assumed to he markovian and satisfy the condition ^(to) = ^o- 
At the same time the function i\A . 1\ ) gives their exhaustive description: 



(A.3) 



where P*-"-* is a density of the probability that the trajectory ^(t) would pass through 
the sequence of intervals ^i, + d^i , . . . , ^n, + d^n at the subsequent moments 
of time ti < . . . < tn, respectively. 

Under these assumptions we can obtain the following representation for an av- 
eraging procedure 



'exp |- £ V, (e(r), eW) dr - V2 {m)]) = J d^e-''^^^'') Q (f, t) , (A.4) 
where d^ = dC,i---dC,n, CLnd the function Q (^^,^',t^ is a solution of the proble 



lem 



OQ 
dt 



t^ta 



(A.5) 



0, 



where \ \ ■ \ \ is a norm in R^. 

Proof. The proof is performed formally under the assumption that all the manip- 
ulations are legal. Denote the left side of the equality ( A.4 ) by / and expand the 
averaging quantity into the Tailor series: 



71=0 



(J-nit), 



(A.6) 



where 



l^n{t) 



Vi{T)dT + V2{t) 



to 



E 



n! 



'q m\{n — m)\ 



vr'^it) { J v,{T)dT 

. to 



(A.7) 
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E 

m=0 



,n — m 



I t Tm T2 ^ 



to to to 

For brevity sake in ( |A.7| ) there was introduced the designation Vi(r) = Vi (^^(r), ^(t) 
V2{t) = V2 {^{t)^ ■ Using the Fubini theorem, we can represent the averaging proce- 
dure in ( p.l8|) as an integration with the weig ht P(") from (^) 

W / dr^ J dr^-i ...J dr^rj . . . V^{n) \ = 

to to to ' 



t T2 

J dij du ... j dii j dTm... j dnP^^^ (el t\u, tjj p(') (4, tj\i 



to 



to 



Changing, where it is necessary, the order of integration we can obtain the following 
representation for the n-th moment ^n{t) '. 



Q (n — m)\ 



d^vr"" (e) Qm (e, e', t) 



(A.8) 



where the countable set of functions Qm is determined from the recurrence 

relations 

t 

Qm (e, e', t)=JdTjdff (e) P^^^ (e', t\v, r) V, {v, Qm-i {v, r) , (A.9) 



m = 0,l,2,.... 



where 



Qo {^ i', t) = p(') (e, t\^o, to) (A.io) 

i.e. in fact the function Qq is independent of Upon the substitution of ( |A.^ ) 



into (|A.6|) we insert the summation procedure under the integration sign and then, 
changing the order of double summation, get the expression 



(A.ll) 
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where 



(A.12) 



n=0 



The representation ( [A. 4 ) is thus obtained. 

It remains to prove that the function Q from ( |A.11| ) is a solution of the problem 
(|A.5| ). Using ( |A.12| ) and ( |A.9D we can easily show that Q satisfies the integral 
equation 

t 

Q (e, i', t) + J drjdfj P(2) (^-^ t\ff, t) V, {ff, i') Q {ff, i', t) = Qo (e, t) (A.13) 

to 

Taking into account that Qo satisfies ( |A.2| ) and the initial condition Qo ('Ci ^o) = 
6 (^^— ^oj and is an integrable function, it can be deduced from ( |A.13| ) that Q 
coincides with the solution of the problem ( |A.5| ). The representation ( |A.4| ), ( [A .51 ) is 
thus obtained. A 
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